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Abstract 

A consistent method for calculating the interquark potential generated by the rela- 
tivistic string with massive ends is proposed. In this approach the interquark potential in 
the model of the Nambu-Goto string with point -like masses at its ends is calculated. At 
first the calculation is done in the one-loop approximation and then the variational esti- 
mation is performed. The quark mass correction results in decreasing the critical distance 
(deconfinement radius). When quark mass decreases the critical distance also decreases. 
For obtaining a finite result under summation over eigenfrequencies of the Nambu-Goto 
string with massive ends a suitable mode-by-mode subtraction is proposed. This renor- 
malization procedure proves to be completely unique. In the framework of the developed 
approach the one-loop interquark potential in the model of the relativistic string with 
rigidity is also calculated. 

PACS number(s): 11.17.+y, 12.38.Aw, 12.40.Aa 

1 Introduction 

Investigation of the quark interaction at large distances is outside the QCD perturbation theory. 
Usually, in this field the lattice simulations and string models are used. 

Calculation of the quark interaction in the framework of string models has rather long 
history (see, for example, papers [ll] - and references therein). In all these investigations 
without exception, only the static interquark potential has been considered. It implies that the 
quarks are assumed to be infinitely heavy. Obviously, this potential, by definition, does not 
depend on the quark mass. Assumption about infinitely heavy quarks is rather crude at least 
for u and d quarks with masses about 140 MeV that is significantly less than the characteristic 
hadronic mass scale ~ 1 GeV. It is clear, that in a general case the interquark potential should 
depend on quark masses. Both the general approach to this problem in the framework of 
QCD []TU[ and the numerical calculations of the light and heavy meson spectra in potential 
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models [[TTl- |]T3| testify to this. Certainly, in this case one should say not about the static 



potential generated, for example, by a relativistic string connecting quarks but simply about 
the interaction potential between quarks having finite mass rather than infinite one. 

The aim of the present paper is an attempt to extend the standard approach to calculation 
of the interquark potential in the framework of the string models |Jl[|, ^ to the case of 
the finite quark mass. It turns out that this program can be realized. To this end, the bound- 
ary conditions in the string model in question should be modified and a new renormalization 
procedure should be developed. 

In a proposed approach, a correction to the string potential due to the finite quark masses 
is calculated both in the Nambu-Goto string model and in the Polyakov-Kleinert rigid string 
model. 

In the Nambu-Goto string with massive ends the quark potential is calculated first in 
one-loop approximation of perturbation theory for arbitrary dimension of the space-time D 
and then via variational estimation in the limit — 2) — ^ oo. As knows, the static quark 
potential generated by Nambu-Goto string in one-loop approximation is compiled by two 
terms, linearly rising confinement potential and universal Liischer l/i?-term |T^. When the 
ends of the Nambu-Goto string are loaded by point-like masses (quarks) then the potential, in 
addition to terms mentioned above, acquires correction dependent on quark mass m and the 
distance between quarks, R. In fact, R and m are involved in potential correction in the form 
of a dimensionless parameter q = (Mo/m)Mo R, where Mq is the string tension. As a result, in 
the correction obtained the limit of small quark mass (Mo/m —>■ oo) is equivalent to the large 
R limit {MqR oo), and the limit of heavy quarks (Mo/m — » 0) is the same as the small R 
limit. At large R (or at small m) the correction to the quark potential reduces to the constant 
determined by m and Mq plus terms of order as R~^ and higher. Thus, in the framework of 
the perturbative calculation, the universal Liischer 1 / i?-term is preserved in the Nambu-Goto 
string with massive ends too. At small R (or in the case of heavy quarks at string ends) the 
mass quark correction to the potential, as one would expect, vanishes. In our calculation, this 
is a direct consequence of the renormalization procedure. It is important that the subtraction 
procedure used for this aim proves to be unique. 

Variational calculation of the potential generated by the Nambu-Goto string with massive 
ends in the limit (D — 2 ) — oo results in the radical expression (see Eq. (4.10)). Quark 
mass contribution to it can be interpreted as a substitution of the critical distance R^ = 
tt{D — 2)/12 Mq from the Nambu-Goto string with fixed ends by quantity R^{m, R) dependent 
on quark mass m and distance R. As a results the radical formula (4.10) has sense not only at 
R > Rc, as in the case of the Nambu-Goto string with fixed ends, but also at R^^^ < R < R, 
where -R^"'^-^ is an effective critical distance determined by condition VlR"^-^-^) = 0. Thus, R^-^-^ 
turns out to be dependent on quark mass m and when m —>■ 0, R^Jf{m) decreases. 

In the rigid string model with massive ends the interquark potential is calculated in one- 
loop approximation. When confining to the quadratic approximation in the string action in 
this model, the dynamical variables (string position vector) can be presented as a sum of two 
terms, u(t, r) = Ui(t, r) + U2(t, r), where Ui(t, r) is a solution to the Nambu-Goto string with 
massive ends and U2(t, r) is additional variable caused by the curvature in the Polyakov-Kleinert 
action. It is remarkable that quark masses only affect on ui(t, r). This essentially simplifies the 
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problem under consideration and enables us to use directly the results for potential derived in 
the Nambu-Goto string with massive ends. In one-loop approximation, the variables ui(t, r) 
and U2(t,r) give additive contribution to the interquark potential generated by rigid string. It 
is true both in the case of the fixed string ends and for the rigid string with massive ends. As a 
result, the mass quark correction to the one-loop potential generated by rigid string is reduced 
to the modification of the contribution from the variable Ui(t, r): the one-loop potential in the 
Nambu-Goto string with massive ends calculated before should be used here. 

In all these calculations we are dealing with an infinite sum of all eigenfrequencies of the 
Nambu-Goto string with massive ends. It substitutes the well know sum over all integer 
frequencies in the Nambu-Goto string with fixed or free ends: (1/2) X^^i ^ = —1/24. For 
obtaining a finite value of this new sum, a mode-by-mode subtraction procedure is proposed: 
each initial eigenfrequency is subtracted by the corresponding frequency of the same string with 
massive ends taking the limit R ^ 0. This prescription turns out to be completely unique. By 
making use of the argument principle theorem, we present this regularized sum in a form of 
the integral. Numerical calculation of the subtracted sum over string eigenfrequencies directly 
and using for this purpose the integral representation give the same result. 

The layout of the paper is as follows. In Section 2 the quadratic approximation for the 
Nambu-Goto string model with massive ends is developed. Upon linearization of the equations 
of motion and boundary conditions, the general solution to them is obtained. The eigenfre- 
quencies of the string oscillations are determined by a transcendental equation. Then the 
canonical quantization is shortly outlined. In Section 3, the interquark potential generated 
by the Nambu-Goto string with massive ends is calculated in one-loop approximation of the 
perturbation theory. In order to remove the divergence, a new subtraction procedure was pro- 
posed. In Section 4, interquark potential, generated by the Nambu-Goto string with massive 
ends, is calculated by making use of a variational estimation of the functional integral in the 
limit when {D — 2) oo. In Section 5, the rigid string model with massive ends is treated. By 
making use of a quadratic approximation for the Polyakov-Kleinert action, the linear equations 
of motion and boundary conditions are derived. Then canonical quantization of this model is 
developed. And finally, interquark potential generated in this string model is calculated in 
one-loop approximation. In Conclusion (Section 6) the obtained results are shortly discussed 
and possible extension of them are proposed. Some mathematical details of calculation are 
presented in Appendices A and B. 



2 Nambu-Goto string with massive ends 

The action of the Nambu-Goto string with point-like masses attached to its ends is written in 



the following way |1I5 
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S = -Ml jj rfS - Y.maj dSa, (2.1) 



a=l 



where ciS is infinitesimal area of the string world surface, Ca [a = 1, 2) are the world trajec- 
tories of the string massive ends and Mq is the string tension with the dimension of the mass 
squared {h = c = 1). 
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For our calculation it will be convenient to use the Gauss parametrization of the string 
world surface: 



a;'^(0 = (t,r;,Ti(t,r),...,a;^-2(t,r)) = 

-(e;u(e)), i = 0,l. (2.2) 

The vector field u^(t,r), j = 1, ...,D — 2 corresponds to D — 2 transverse components of x^, 
while the ^ = (t,r) are the coordinates on the string world sheet. The infinitesimal area is 
given by ciS = dt dr y^, where g is the determinant of the induced metric on the world surface 
of the string, g^j = diX'^djX^. The metric of the D-dimensional space-time has the signature 
(+,-,...,-). 

In this parametrization, the induced metric gij in the quadratic approximation has the 
following components 

9ij = ^ij - UiUj , (2.3) 
where uu = J2f=i'^'''^J'''- From Eq. (2.3) we obtain 

g'^ = g-'[il-nl)S,j + u,u,-] , (2.4) 

g = det(^i,)~l - (2.5) 
(here Uq = du/dt — u, vii — du/dr — u'). The hne elements dsa, a — 1,2, take the form: 

dsa^[l-\ix\t,Ta)]dt. (2.6) 

After neglecting unimportant constants, the action (2.1) becomes 

S ~ -Ml {t-2 - ti)R + Jdtjdr [u\t, r) - u'\t, r)\ + 

ti 

2 *2 

+ E ^/ dtu\t,ra), ri = 0, r2^R. (2.7) 



tl 



The last two terms in Eq. (2.7), important for deriving dynamical equations, can be rewritten 
as follows 



S c:, '-^ J dt J dr [u'it, r)e{r) - u\t, r)] , (2.8) 



tl 

where e{r) is the mass density distribution along the string 



sir) = 1 + ^ m + d{R - r)] . (2.9) 



For semplicity we assume that the string ends are loaded with equal masses, mi = m2 — m. 
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Equations of motion and boundary conditions can be deduced from (2.7) or (2.8) by relations 

d ( dLstr\ d ( dLstr\ „ 

and 

^--1 7^7r^ = 0' «=1'2, r, = 0,i?. 2.11 

Lstr is the string Lagrangian density and L^, a = 1, 2 are the Lagrangians of the massive ends 
of the string in Eq. (2.7). Equations of motion are given by 

□u = 0, (2.12) 

where □ = d'^/dt^ — d'^/dr'^, while for boundary conditions we find 

mil = Mgu', r = 0, (2.13) 

mii = -M^u', r = R. (2.14) 
Equation (2.12) admits the following solutions: 

u^{t, r) ~ exp[i{iu / R)t] u{r) , j = 1, 2, . . . , D - 2 , (2.15) 

in which the string length R has been introduced in order that uj be dimensionless. By substi- 
tuting (2.15) into (2.12), (2.13) and (2.14) one obtains 

u"(r) + ^u(r) = 0, (2.16) 

uj^u{0) = -qRu'{0) , (2.17) 
uj'^u{R) = qRu\R) , (2.18) 



where g is a dimensionless parameter 



,= M^. (2.19, 



m 



The differential equations (2.16)-(2.18) are linear, therefore each component of the transverse 
oscillation of the string satisfies the same equations. Hence, we can write the general solution 
as a superposition of plane wave solutions 



\/2 

u'{t, r)=tj-J2 exp[-t{iUn/R)t] Un{r) , J = 1, 2 , (2.20) 



where the amplitudes (Fourier coefficients) obey the usual rule of complex conjugation, 
«n = «!„,, in order that be real. The eigenfunctions Un{r) in (2.20) are defined by 



uJr) = N„ 









cos (u^^) 


— — sm 







(2.21) 
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where Nn are normalization constants and the eigenfrequencies Un are the roots of the tran- 
scendental equation 

tanu; = — -. (2.22) 
ijj^ — 

On the C(j-axis these roots are placed symmetrically around zero. Hence they can be num- 
bered in the following way cuq = 0, uj-n = —uJn, n = 1,2,.... Therefore it will be sufficient 
to consider only the positive roots. The eigenfunction u„(r) obey the orthogonality conditions 
with the weight function e{r) |17| 

R 

drun{r) Um{r)e{r) = R6nm, (2.23) 



while the functions u'^{r) satisfy the usual orthogonality conditions 



dru^{r)u^{r) = ^S^m, (2.24) 



where the eigenfrequencies Un are solutions of Eq. (2.22). 

The density of the canonical momentum 'p'{t,r) is defined in standard way 

P'{t,r) = -^ = Mlu^{t,T)e{T), (2.25) 

in which Ltot is the Lagrangian density in action (2.8). From (2.25) the total momentum of the 
string is given by 

R 

pi it) = J drjP{t,r). (2.26) 



Obviously, in the problem under consideration, we can put Pj = 0. The canonical Hamiltonian 
is defined by 



R 

H = jdr[p{t,r)\i{t,r) - Ltot\ 



M2 



= ^Jdr [ii'it, r)e{r) + n\t, r)] , (2.27) 



In terms of the amplitudes anj it reads 

1 oo D-2 

^ = R E E <^ + <^ <) ■ (2-28) 



n=l j=l 

When we quantize, ui{t,r) and its conjugate momentum p^{t,r) become operators with 
canonical commutation relations 

[u\t,r),jP{t,r')] =i5'H{r-r'), (2.29) 
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This implies that the Fourier coefficients become operators and satisfy the relations 

[<,<] = u;„(5*^(5„+^,o, (2.30) 

i, j = 1, . . . , D — 2 , n,m = ±1, ±2, .... 

The creation and annihilation operators, respectively and a„, in Fock space are intro- 
duced in usual way 

(2.31) 
(2.32) 



at 



j + 



nm 1 ^5 ^ 1)2, 



and through them, the Hamiltonian (2.28) takes the form 

oo D—2 2 1 °° 

^ n=l j=l 2 i< 

The last term in (2.33) is the usual Casimir energy [|l^, ^ 



(2.33) 



3 One— loop potential generated by Nambu— Goto string 
with massive ends 

In this section we shall investigate of the interquark potential arising from the action (2.1), 
via the perturbation calculation. For this purpose the Euclidean version of the model under 
consideration should be treated. 

Up to the second order in u the string action (2.1) in Euclidean space is given by (ti = 
0, t2 = T) 

" 2 " 



T R 

= fdt fdr 







1 + -ud'^u 
2 



(3.1) 



ri = 0, r2 = R, 

where = d'^/dt'^ + 8"^ /dr"^. It generates one-loop Feynman diagrams in perturbation theory. 
The potential V{R) between massive quark separated by a distance R is defined in terms of 
functional integral as follows [pO| , pl[| : 



-TV(R) 



[Dvl] e 



-5b[u] 



T ^ OO . 



(3.2) 



From action (3.1), after functional integration, the interquark potential reads 

1 



ViR) = lim _ 



^ ^ D-2 
dt drM^ + TrlnG"^ 



.0 



(3.3) 
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In Eq. (3.3) G ^ = 9^ is the operator generated by the quadratic part of the action Se 
and G is the Euchdean Green function 

d^G{x,x') = 6{x,x'). (3.4) 

In the momentum space this function is written as 

/to "T fi^n 

where fc„ are admissible values of the wave vector for the field u(t,r), determined by the 
boundary conditions (2.13) and (2.14) 

kn = ^. (3.6) 

As before Un are the roots of frequency equation (2.22). By using the definition of the Tr in 
the momentum space 

auj 



Tr... = T/— , u; = ko, (3.7) 

where the sum is spread over all the discrete values of the component kn, we can present the 
functional trace in Eq. (3.3) in the following way 

■ du! 



IYlnG-^ = rE J^Hu^ + kl). (3. 



In the analytical regularization the cj-integration can be performed by the formula P] 

fpHcu' + a') = V^. (3.9) 



It follows that up to the one-loop level the interquark potential (3.3) is given by 

V{R) = MlR+ ^^v{q), (3.10) 

in which 

^(9) = ^ E ^n{q), (3.11) 

^ n=l 

and ujn{q) are the roots of Eq. (2.22). 

The sum in (3.11) diverges |19] so that it is necessary to develop a subtraction procedure 
physically acceptable in order to get the renormalized function f''^"(g) in a unique way. We 
propose here a mode-by-mode subtraction procedure. To this end we consider the asymptotic 
of eigenfrequencies w„(g) for large n but fixed q 



w„(g) ~ nTT + ^ - fl + I) M! + o(n-5), n^oo. (3.12) 
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The divergence of the sum (3.11) is due to the first two terms in (3.12). At the same time we 
consider the asymptotic forms of Eq. (2.22) at — > 



and at 5 — > oo 



tana; = ^ , (3.13) 

UJ 



cot a; = (3.14) 



The asymptotic of the roots of these equations for large n and fixed q are given, respectively, 

by 



\2 



u{q) ^n7r+;^-fl + ^) M_ + o(^-5) , ^ ^ oo , (3.15) 

and 



Thus we see that in order to remove the divergencies in sum (3.11) we must substitute a;„ by 
ujn — oOn- The roots of Eq. (3.14) cannot be used for this purpose in view of their asymptotic 
(3.16). Therefore the subtraction procedure proposed here is unique. 

In terms of the function v{q) the proposed renormalization procedure is defined in the 
following way 

^""(?) = v{q) -Ii^,_.o^(?) + ^^0. (3.17) 



The sjTnbol \\va.q^Qv{q) means that we must take the asymptotic expression of the function 
v{q). Analogous renormalization scheme is used by calculating the Casimir energy in some field 
models with special boundary conditions. We have introduced here a constant Vq in order to 
satisfy the requirement that aX m = oo Eq. (3.17) must give v^^^^^q) for string with fixed ends, 
i.e. 

v''^{q) = v{q) - mv{q) + v{q)\m = oo ■ (3.18) 
The frequencies of the string with fixed ends are 

cUn — mr . (3.19) 
Therefore the last terms in (3.18) is defined by 

TT °° 

Vq = v{q)\m = oo = Z] ^- (3-20) 

^ n=l 

This sum can be renormalized by making use of the Riemann zeta-function ((s) = J2^=i n~^: 

oo 

v{q)\m = oo = 2 = 2^^"^^ ^ "24- ^^"^^^ 

The final expression for the renormalized function v^^"'{q) is given by 

TT . 1 

n=l 



^""(^) = - ^ + o E Mq) - c^niq)] . (3.22) 
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This function has been calculated numerically. As shown in Fig. 1, v^'^^{q) is a monotone 
increasing function of the dimensionless parameter q. 

The interquark potential (3.10) in terms of the function reads 

V{R) = MlR+ ^^—^ v'^'^iq) . (3.23) 
R 

Equation (3.23) shows that the potential generated by string with massive ends depends on 
the quarks masses through the dimensionless parameter q = M^R/m. Taking into account 
this fact we can give here a very clear physical meaning of the subtraction procedure proposed 
above for obtaining renormalized quantum correction to the potential (3.22). The limit g — *^ 
may be obtained by putting i? — 0. Therefore, v^^^{q) is determined as a difference of v{q) in 
two points: r = R and r — > 0. In the same way, the potential (3.23) is equal to the work that 
should be done for removing quarks connected by string to the distance R, provided that at 
the beginning they were very closed to each other. 

Certainly, this renormalization procedure should be followed by substitution for Mq, the 
renormalized (physical) value of the string tension. For simplicity we do not introduce new 
notations and will keep this point in mind. 

Now we propose a different way for calculating the sum in Eq. (3.22). It utilize an analytical 
method that employes the following integral formula from the complex analysis |^ . 



Let us consider an analytical function f(z) with zeroes of order Uk at points z = ak and 
with poles of order pi at points z = bi in a. region bounded by a contour C. From Cauchy's 
theorem it follows that 

-L IdzzQ^ = ^ hzz[\n fiz)]' = n,a, - ^ ph. (3.24) 

In order to obtain the renormalized function v'^^^{q) through Eq. (3.24), one has to choose 
an appropriate expression for the function f\z). Equation (3.18) implies that f{z) should be 
taken in the following form 

f{uj) = smu. 3.25 

— Iquj cot uj 

Here the numerator is in fact Eq. (2.22) and the denominator is the same equation in the limit 
g — s> (Eq. (3.13)). The multiplier sincu is introduced to take into account the frequencies 
(3.19). 

When applying the formula (3.24) to nominator, denominator and multiplier sincj in (3.25) 
separately with the counter C encircling the real positive semiaxis in the complex plane cj, the 
contribution of the poles into (3.24) is obviously absent. Now we deform the contour C so that 
it transforms into the imaginary axis on the complex plane. After integrating by parts, the 
final expression for the renormalized function v^^'^{q) becomes (Appendix A) 

v'^^'il) = - J + ^(^) ' (3-26) 
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where 

1 



/(g) = ^ Jdy In 



^ (v/q)'^ + 2(y/g) coth y 



(3.27) 



The integral /(g) has been calculated numerically for different values of the parameter g. The 
results obtained for f'"'^"'(g) by making use of (3.22) and (3.26), (3.27) are certainly the same, 
(see Fig. 1). 

Finally we analyze the asymptotic expansion of the potential V{R) in the limit g — * oo or 
R —>■ oo. In this region the renormalized function v^'^'^{q), Eq. (3.26), is given by (see Appendix 
B) 

v-{g)^-f^ + —g+^, g = ^, g-oo. (3.28) 
24 TT 27rg m 

At large R the interquark potential (3.23) can be written up to second order in 1/R as follows 

, (/} - 2) ln2 M2 7i(D - 2) 1 

The first terms in (3.29) is the usual linearly rising term Mq R (confining potential). The 
second term is a constant dependent on the quark mass. The third term is the universal Liischer 
l//?-term and the last one is a correction to the string potential due to the finite quark mass. 



4 Variational estimation of the string potential 

In preceding Section the string potential has been calculated by using the perturbative theory 
for an arbitrary dimension of space-time D. Otherwise this potential can be derived in the 
limit D ^ oo hj making use of the variational estimation of the functional integral 
Let us turn to the initial equation (3.2) determining the string potential 

e-^^(^) = J [Du] e-^-M , T oo , (4.1) 

where Se is the Euclidean action 




The l/(/} — 2)-expansion is carried out in standard way [l|]. Let us introduce the composite field 
aij for diudjU and constrain aij = diudju through the Lagrange multiplier a*-'. By using the 
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exponential parametrization of the 5-function, with the understanding that the a*-' functional 
integrals run from —ioo to +ioo, and calculating the u-functional integral, Eq. (4.1) becomes 



(4.3) 



where 



Ml jdtjdr 





det {6ij + aij) - -a'^ ai 



+ 



(4.4) 

The boundary terms in (4.2) should be taken into account under finding the eigenvalues of the 
differential operator —diO^^dj in (4.4). The 1/{D — 2)-expansion is generated by expanding 
the action (4.4) around its stationary point. We shall assume that at the stationary point a^^ 
and (Jij are diagonal matrices, a and a respectively, independent on t and r []T]]. From Eqs. 
(3.6), (3.7) and (3.9) we obtain 



Trln(-a,a^^a,) =T^ / 



11 7.2 ^ 



2T a 
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v{q) 



R V 

where v{q) is defined by (3.11). Now, action (4.4) can be rewritten as 



MlRT 



1 



Vl + (JoVl + o-i - 2 ('^°'^o + aVi) 



+ 



(4.5) 



(4.6) 



with a" = a*, an = ai, i = 0, 1. Equating to zero the variations of the action (4.6) with 
respect to ctq, ai, and solving the corresponding set of algebraic equations allow us to 

find the stationary points 



a 



Vl - 2X, 
X 



1 



a 



1 - 2A ' 

where the dimensionless function A(g) is defined by 



c^i 



= -A, 



(4.7) 
(4.8) 

(4.9) 
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Here we have used the renormahzed function v{q) from Eq. (3.22). 

Putting (4.7) and (4.8) into (4.6) yields the string potential to leading order in 1/{D — 2): 



V{R) = M'.rJ 1 - (l - ^-^I{q)y 



(4.10) 



where — 7r{D — 2)/12Mq and I{q) is defined in Eq. (3.27). The quark mass contribution 
in this formula can be interpreted as a substitution of R^ from the Nambu-Goto string with 
fixed ends by function R?'{m, R) dependent on quark mass m and distance R: 

24 



R\m,R) = Ri 



TT 



'MlR? 
m . 



(4.11) 



This results in an extension of the applicability of Eq. (4.10) to i? < Re- More precisely, 
formula (4.10) has sense not only for R > Rc, that takes place for the string with fixed ends, 
but also in the region 

Rl^^ < R < Rc, (4.12) 

where R^-^^^ is defined by equation 

R'{m,Rlf^) = {Rlf^f. (4.13) 

At -R = R^^'^ the interquark potential (4.10) vanishes. Now critical distance R^^^ is determined 
by the quark mass m and when m — R'^^^{m) also decreases. In Fig. 2 the dimensionless 
potential V{R)/Mq as a function of the dimensionless distance p = Mq R is plotted for different 
values of the ratio = Mo/m. 

At the end of this Section we present an asymptotic expression for the interquark potential 
(4.10) at -R ^ oo. The asymptotic expression for function I{q) when q — > oo, (or R — > oo), is 
given by (B.5) 

, ln2 M^R 1 m 
7(g) ~ ^ + — 



TT m 

Taking into account (4.14), we obtain 



27r R ' 



R 



oo . 



(4.14) 



V{R) r^M^R + 



{D - 2) In 2 Tx{D - 2) 



TT 



+ 



D-2 m 
27ri?2 



1 + 



m 24: R 

7r{D-2) In 2 



1 + 



12 {D - 2) (In 2)2 M2 



+ 0{R-') , 



R 



+ 



oo . 



(4.15) 
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The interquark potential calculated via variational approach, at large values of R has some 
additional terms as compared with the potential (3.29), obtained by perturbative method. In 
particular, it turns out that the Liischer l/i?-term 

tt(D - 2) 



24 R 



should be now substituted by 

7r{D 



24 R 



1 + 



12 (D - 2) (In 2)2 M2 



(4.16) 



(4.17) 
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5 Quark mass corrections to the rigid string potential 



In this Section we calculate the mass quark corrections to the one-loop interquark potential in 
the framework at the rigid string model As known, this model can be treated as an 

effective one taking into account the finite thickness of gluonic tube [^-p^. The basic aim of 
this calculation is to show the principal applicability of the proposed method to the rigid string 
model with massive ends. Variational estimation of the interquark potential in the framework 
of this model will be published elsewhere. 

The Polyakov-Kleinert action for rigid string with massive ends has the form 



S = -Ml 



a 



1 - -r^Ax^Ax^ 



2 r 
^rUa I dSa 
a=l i 



(5.1) 



where the new parameters and a are, respectively, the radius of gluonic tube and a dimen- 
sionless constant; A is the Laplace- Beltrami operator for the induced metric gij 



A 



1 



d 



-9 



d_ 



(5.2) 



In the Gauss parametrization (2.2), the operator (5.2), up to the second order in u, can be 
written as 

A ~ □ + 0{v?) , (5.3) 

Now action (5.1) reads 



S 



t2 R 

-Ml fdt [dr 



2 2 2 ^ 



2 *2 

- Jdt^\t,ra), n=0, r2 = R. 



a=l 

The equations of motion and boundary conditions for the action (5.4) are 

(1 + ar^D) Du = 0, 
m 



1 + « □) u' = ^ ii , r 

m 

^ ii' = — 



:i + ar?n)u' = - — 2 u 
Ml 



0, 
r = R. 



□ u 



0, 



r 



0,R, 



(5.4) 

(5.5) 
(5.6) 

(5.7) 
(5.8) 



(mi = 777-2 = ITT')- The Lagrangian in the action (5.4) depends on the first and the second 
derivatives of the string coordinates, therefore the number of obtained boundary conditions is 
twice compared with the Nambu-Goto string. 
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The boundary value problem (5.5)-(5.8) reduces to the two independent boundary problems. 
Indeed equations of motion are given by product of commuting differential operators (1 + a □) 
and □. Hence, the general solution to this equations can be represented as a sum of two terms 



u(t,r) = ui(t,r) + U2(t,r) 



(5.9) 



where 



□ ui = 

m 

m 



2 ill , r = , 



u 



i-+^ui, r = R. 



and 



1 + a<n)u2 







(5.10) 
(5.11) 

(5.12) 

(5.13) 
(5.14) 



U2(t,0) = U2{t,R) = 0. 

In this case, ui(t, r) is the solution for the Nambu-Goto string with massive ends that we 
have analyzed in Section 2. The string rigidity is taken into account by function U2(t,r). The 
general solution to Eq. (5.13) obeying (5.14) can be presented as 



V2 ^ 

* -TV y 6Xp 



0^ 



Vn[r) 



J = 1,2, 



D 



(5.15) 



The eigenfunctions Vn{r) are given by 

Vn{r) = -V-n 



N„ sin nn — 
R 



n 



1,2, 



(5.16) 



where N'^ are normalization constants. For the natural frequencies, in (5.15) we have 



— 1/_ 



(7™)2 + 



art 



n 



1,2,. 



(5.17) 



The amplitudes satisfy the usual relations of complex conjugation /?* = /?_„, n = 1, 2, . . . . 

The Hamiltonian formulation of the model under consideration is developed in the following 
way. According to Ostrogradskii PBI, EU[ the canonical variables are defined by 



q{ = 



Pi 



dL 



tot i i 



where Ltot is the Lagrangian density in action (5.4) 



qi = u\ 



^^T, J = 1,2, 



D 



(5.18) 
(5.19) 



Mi 



^tot 



eir) u 



u'2 - ar?(nu)^ 



(5.20) 
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Putting Ltot and (5.9) into (5.18), (5.19) and taking into account Eqs. (5.10) and (5.13) one 
obtains 

qi = Ui + U2 , q2 = ui + U2 , (5.21) 
pi = [e{r) + arln]h, pa = - a M ^ □ u = U2 . (5.22) 
The canonical Hamiltonian is defined by 

H = Jdr [piqi + paqa - Lt^t] , (5.23) 



In terms of Fourier amplitudes it becomes 

1 oo D-2 1 oo D-2 

H = TT^Ei: («n,<- + <•«„,) - iPn.P:, + (3t,Pn,) ■ (5.24) 

n=l j=l n=l j=l 

The quantum theory is based on the canonical commutation relations 

Kit,r),pi{t,r')] = tS,,5'^ Sir-r'), a=l,2, t, j = 1,2, . . . , D - 2 , (5.25) 
or in terms of the Fourier amplitudes 

[a^,a^] = Un S'^ 6n+m,o , [Pn,PL] = i^n ^n+mfi , n,m = ±1,±2,... . (5.26) 
By introducing in standard way the annihilation and creation operators 

< = (c.„)-i/^<, = (c.„)-^/2a^, 

K = (l^nr'^'Pi, V: = (5.27) 

n = 1,2,..., 1,3, = l,2,...,D-2, 
the Hamiltonian operator (5.24) acquires the following form 

1 oo D-2 1 oo D-2 

^ = o E E <^< - o E E + 

^ n=l j=l ^ n=l j=l 

T-) _ 2 / oo OO \ 

+ ( E - E \ ■ (5.28) 



The last two terms in (5.28) define the Casimir energy in the model under consideration [RD 



It is important to note that the second oscillation mode with frequencies i/^, responsible for the 
string rigidity, gives a negative contribution to the energy as compared with the oscillation of 
the first mode with frequencies uo^. This is also true for the Casimir energy (see the last two 
terms in (5.28)). It is a direct consequence of the classical expression for the total energy in the 
rigid string model (5.24). We point out that this defect is typical in all the theories with higher 
derivatives. To remove it, certainly in formal way only, the quantum states with negative norm 
are used sometimes |53[ . 



16 



Now we calculate the rigid string potential in one-loop approximation. Again we shall treat 
the Euclidean version of the model under consideration. 

The interquark potential is given by Eq. (3.2) with the Euclidean action {ti — , t2 — T) 



T R 

Ml jdt jdr 





1 + - u(l + arid')d^n 



(5.29) 



As in section 4, the boundary terms in action (5.4) will be taken into account by finding the 
proper eigenvalues of the corresponding differential operator. 
After functional integration, the potential takes the form 



V{K) 



lim — 

T-*oo T 



T R 



jdt jdr Ml + ^^^y^TrlnC"^ 



.0 



(5.30) 



where the operator G ^ = (1 + a 9^) 5^ corresponds to the inverse of the Euclidean Green 
function 

(1 + a rl d^) G{x, x') = 5{x, x') . (5.31) 
In view of (3.7), the functional trace in Eq. (5.30) can be written as 



TrlnG"^ = Trln^^ + Trln(l + arld"^) = 
TEjZ + k'n'')+Ti: j^ In [1 + arl {u^ + k^r- 



(5.32) 



k[}^ are the admissable values of the wave vector for the field Ui(t, r) determined by the bound- 
ary conditions (5.11) and (5.12), fc^^^ = Un/R, n = 1,2,..., and k!^^ are the admissable 
values of the wave vector for the field U2(t,r) determined by boundary conditions (5.14), 

k^^ — Vn/R, IT- — 1,2, By using the analytical regularization, formula (3.9), Eq. (5.32) 

can be reduced to 

OT^ rp OO 

TrlnG-^ = —v{q) + -7i + e"! , (5.33) 



R 



n=l 



where the dimensionless parameter e is defined by 



(5.34) 



and the function v{q) is given in (3.11). Now the interquark potential (5.30) assumes the form 

D-2 



V{R) = MlR + 



R 



^ n=l 



(5.35) 



In Eq. (5.35) the second mode of oscillation gives a positive contribution to the energy. As was 
mentioned above, it means that the formalism, applied here, effectively uses for excitation of 
the second mode quantum states with negative norm. 
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To remove the divergences in Eq. (5.35) one has to use a subtraction procedure. The function 
v{q) has been renormahzed in Section 3 (Eqs. (3.22) and (3.26)). Now we renormahze only the 
last sum in (5.35). It can be rewritten in the following way 



1 oo 1 / -i-oo 1 \ 

w{e) = - E + e-i = t E Vn^ + e-i - ^ , (5.36) 

^ n=l ^ \n=-oo Ve/ 

where the last term removes the term with n = in the sum. The renormalized function 
w'^^"'{e) is obtained from equation (5.36) by subtracting its value for infinite string 

w'''''{e) = w{e) - lim w{e) = 
Equation (5.37) can be expressed in terms of the function 



S{x) = - E - /^^ +x = E -^i(27™v^), (5.38 

^ \n=-oo / ^ n=l 



where Ki(z) is the modified Bessel function ||3J]. This representation for S{x) is convenient for 



investigating its behaviour at large x, because the modified Bessel function decays exponentially 



3^ . For small x another representation of the function S{x) was proposed By using the 



definition of the function S{x), Eq. (5.37) becomes 



«;-(e) = ^Sie) - (5.39) 



Finally one-loop calculation of the interquark potential gives 



V{R) = M^R + — 



R 



,ren 



IT 1 ~ 1 _ f27m 



(5.40) 



which depends on the quark mass and string rigidity through the dimensionless parameters q 
and e, respectively. 

In the asymptotic limit, R — * oo, the string potential (5.40) assumes the following form 

V(R) ~ + m - 2 ^ ^ — — ^ =^ + 

^ ' ° ^ ' \ n m 2AR 2V2^rJ 



+ - - 2) ^M-2V2 R/V^r.] , R^oo. (5.41) 

The first term in the right-hand-side of (5.41) is linearly rising potential, MqR; the second 
term is a constant determined by quark mass m, string tension Mq and coupling constant y/a 
in the rigid string model action (5.1); the third term is the universal Liischer l/i?-term, the 
fourth term is the mass quark correction, [D — 2)m/27iMQ R^. The last term in (5.41) vanishes 
exponentially when R — »• oo. For different values of the parameters Mo,m, a and the total 
constant term in (5.41) may, in principle, change sign. 
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6 Conclusion 



In this paper we have developed a consistent method for calculating the interquark potential 
generated by relativistic string with point-like masses (spinless quarks) at its ends. The ob- 
tained results indicate that the correction to the potential due to the finite quark masses turns 
out to be considerable near the critical (or deconfinement) distance. When the quark mass m 
decreases the contribution of this correction rises. However, the formula obtained cannot be 
obviously used at very small m. The point is that all our corrections are based on the linearized 
dynamical equations in the string theory which are, in some sense, equivalent to the nonrela- 



tivistic approximation in the initial string action [[Iq , pO| . At small m, when large velocities of 



the string ends are important]^, nonrelativistic approximation is certainly not applicable. 

The extension of the proposed method for investigating the model of the relativistic string 
with massive ends at finite temperature is of undoubted interest. 
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APPENDIX A: COUNTER INTEGRAL REPRESENTATION FOR THE SUM 
OF THE STRING EIGENFREQUENCIES 



Here we calculate the sum (3.20) by making use of the counter integral (3.24). Putting in 
(3.24) f{z) = sincj and deforming counter C as was explained in Section 3, we obtain 



^ +00 



n=l 



J = - J2 n = -— dyy cothy. (A.l) 







The integral (A.l) diverges at the upper limit. To regularize it we introduce into the 
integrand the cutting multiplier e"*^^, e > 0. After doing the integral we take the limit e — >■ 

jren ^ ^-^ j .^ 2) 

27r ^^0 " ^ ^ 



^It is well known |l5[ that free (massless) ends of the Nambu-Goto string are permanently moving with light 
velocity. 
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where 



ttC 2,- 



and ({z, q) is the Hurwitz zeta-function [g4 



(A.3) 



(A.4) 



^ 1- 

+ lim 

24 ^-0 



1 



^0 {n + qy 

Substituting (A.3) into (A. 2) we obtain 

jren 

The pole in (A. 5) should be dropped as it is usually done in analytical regularization. Hence 

jren = _ _ . (A.6) 



27re2 



(A.5) 



TT 

24 



Application of (3.24) to the fraction in (3.25) alone results, after deforming the counter C 
and integrating by parts, in (3.27) without any divergences. 



APPENDIX B: ASYMPTOTIC EXPANSION OF THE INTEGRAL /(g) 



In this appendix the asymptotic expansion for g ^ oo of the function /(g) entered in f ^^"(g) 
will be obtained. In terms of the variable x = q'^ y, Eq. (3.27) reads 



CO , \ 

Uq) = 1- [dxlnil + — I . (B.l) 

271 J \ x^ + 2xcothqxJ ^ ' 



Let us divide the range of integration into two regions, (0, 1/g] U [1/g, cxd), with g oo. The 
integral in (B.l) can be written as 

. / i \ 

+ 



/(g) = 1- jdxlnll + 

^ 

oo / 



+ 2x cothgx 

- [dx In (1 + — , 

271 J \ x^ + 2x cothgx 

1/9 ^ 



oo 

a 

+ 



Z7T 271 

For X G (0, 1/g] we can substitute cothgx by (gx)^^ and the integral Ji(g) in (B.2) gives 

J,(,)./ll„(l + -^).ll„| + l. (B.3) 
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When X & [1/q, oo) the function cothgfx can be approximate by 1, so that the integral J2{q) in 
(B.2) assumes the following form 

oo 

J,fe)./d.ln(l + -^).21.2--ta|. (B.4) 

V? 

Putting Eqs. (B.3) and (B.4) into (B.2) yields 

, , ln2 1 , , 

/(g) ~ q + - — , g^oo. (B.5) 

TT Ztt q 
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Figure Captions 



Fig. 1. Renormalized sum over the natural frequencies of the string with massive ends, 
v'^'^^{q), calculated by making use of the definition (3.22) and the integral representation (3.26), 
(3.27). The results are, certainly, the same. 

Fig. 2. Dimensionless interquark potential V{R)/Mo calculated through Eq. (4.10) for 
different values of ratio /j = Mo/m. The Alvarez result [1] for the Nambu-Goto string with 
fixed ends is obtained at /i — 0. 
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